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During  the  previous  and  current  funding  periods  by  the  AFOSR  the  following  two  main  areas  I  and 
II  of  research  were  pursued: 

I.  Numerical  solutions  of  strongly  nonlinear  partial  differential 
equations  of  Monge- Ampere  type  and  their  applications. 
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Below  we  describe  progress  made  on  these  topics. 

Part  A:  General  Remarks 

The  physical  phenomena  described  by  nonlinear  partial  differential  equations  have  become  at  present 
the  central  theme  of  investigations  by  many  researchers.  A  good  understanding  of  most  physical  processes 
requires  accounting  for  nonlinear  effects  and,  consequently,  methods  for  studying  nonlinear  equations  have 
to  be  developed. 


Among  nonlinear  equations  the  Dirichlet  problem  for  the  Monge-Ampere  equation  is  the  model  case 
for  fully  nonlinear  equations.  The  problem  is  formulated  as  follows. 

In  Euclidean  plane  R 2  with  Cartesian  coordinates  x,y  consider  a  bounded  domain  Q,  a  nonnegative 
function  / :  Q  -*  [0,<»),  and  a  continuous  function  3D  -» R.  It  is  required  to  investigate  solubility  of  the 
problem 


M(z) 


B2z  B2z  (  d2!  \  ..  n 

dx 2  By2  ~'BxBy  '  ’ 


(D 


*  Ida  =  $ 


(2) 


The  equation  (1)  is  a  representative  of  the  class  of  nonlinear  equations  of  Monge-Ampere  type.  Such 
equations  have  been  studied  by  many  authors,  mainly  in  connection  with  problems  of  existence  and 
uniqueness  of  surfaces  with  prescribed  metric  or  curvature  functions.  However,  they  also  have  other 


important  applications.  In  particular,  the  leading  term  in  the  "balance  equation”  in  dynamic  meteorology 
has  the  form  (1).  In  a  more  complicated  form,  an  equation  of  this  type  appears  in  the  von  Karman  system 
of  equations  for  elasticity  and  also  in  inverse  problems  of  geometric  opdcs.  It  also  turned  out  that  recent 
progress  in  the  study  of  fully  nonlinear  equations  became  possible  after  important  properties  of  the  equa¬ 
tion  (1)  were  discoveredtFor  quasilinear  elliptic  and  parabolic  equations  the  use  of  equations  of  the  form 
(1)  is  crucial  in  obtainUlg  C°  and  Holder  estimates. 

Part  B:  Numerical  solution  of  equations  of  Monge-  Ampere  type 


In  view  of  the  important  practical  applications  several  heuristic  approaches  were  suggested  earlier 
for  numerical  soludon  of  some  modified  forms  of  (1),  (2).  Though  no  rigorous  analysis  of  these  methods 
exists,  one  may  note  that  they  all  are  local  methods  based  on  a  finite  difference  approximadon  and  lineari- 
zarion.  Because  of  the  strong  nonlinearity  of  M  this  approach  might  be  successful  only  in  a  neighborhood 
of  the  true  soludon  and  therefore,  if  a  priori  good  initial  approximadon  is  not  available,  these  methods  will 
not  produce,  generally,  a  sequence  converging  to  the  true  soludon. 

We  have  investigated  this  problem  in  detail  and  obtained  the  following  results: 

•  A  special  discretization  scheme  for  (1),  (2)  was  suggested,  different  from  standard  finite  element  or 
finite  difference  schemes.  It  can  be  shown  that  latter  ones  in  known  forms  will  not  work  here. 


•  an  iterative  method  has  been  developed  for  solving  the  discretized  version  of  (1),  (2); 

•  the  question  of  finding  an  initial  approximation  in  our  scheme  is  completely  and  effectively  resolved;  it 
is  just  a  routine  step  of  the  iteration  process; 

•  the  iterations  are  self-corTecting; 

•  global  convergence  is  established; 

•  our  algorithm  is  suitable  for  a  parallel  computer; 

•  a  computer  code  has  been  written  and  tested  on  a  serial  and  parallel  machines. 


•  a  graphics  package  accompanying  the  code  was  also  developed. 


The  experience  gained  in  testing  our  procedure  on  different  types  of  examples,  including  ones  with 
large  gradients,  shows  that  its  most  effective  use  is  in  combination  with  some  fast  Newton-type  scheme. 
Consequently,  a  particular  criteria  is  established  for  checking  when  the  current  approximation  can  be  used 
as  the  beginning  step  for  a  converging  Newton-type  iterative  procedure. 


•  In  this  combined  scheme  we  proved  convergence  (quadratic)  of  the  Newton  iterates. 


ftp 

I  ~ 


The  computer  code  for  the  method  is  quite  sophisticated;  it  involves,  as  a  step,  construction  of  a  con-  ^ 
vex  hull  of  sets  of  points  in  R3.  There  are  here  different  aproaches  and  the  effectiveness  of  the  algorithm  ( 
depends  on  it  We  have  tested  various  schemes  and  the  code  presently  is  a  substantial  improvement  over  *  __ 
its  original  version  of  1984-85.  Recently  our  work  here  was  followed  by  Baldes  and  Wohlrab  in  Germany. 
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Part  C:  Numerical  solution  of  the  "balance  equation"  arising  in  weather  forecasting 


We  developed  extensions  of  our  numerical  procedure  to  elliptic  equations  of  the  form 


32z 

dx2 
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A  special  case  of  this  equation  is  the  "balance"  equation 

+<a*((-Wgnd«)-/  (4) 

to  be  solved  in  a  convex  bounded  domain  QcR2  subject  to  the  boundary  condition 

zldn  =  <t>.  (5) 

where  <p  is  a  continuous  function. 

This  equation  relates  the  rotational  component  of  the  wind  velocity  and  the  geopotential.  It  is  used  in 
a  model  of  numerical  weather  prediction. 

For  the  problem  (4),  (5)  we  define  a  notion  of  a  generalized  solution  by  replacing  (4),  (5)  by  an 
equation  in  measures 


)W+I^(x)^z)dxdy  =  lfdxdy 

where  to  is  an  arbitrary  subdomain  of  Q  with  piecewise  smooth  boundary.  The  first  term  on  the  left  fits 
into  the  scheme  developed  by  us  earlier.  For  the  second  term  we  construct  a  special  discretization 
compatible  with  the  one  used  for  the  first  term.  As  a  result  we  obtain  a  system  of  quadratic  equations.  This 
system  is  solved  by  an  algorithm  which  is  an  extension  of  the  one  in  part  B.  A  code  for  this  procedure  has 
been  written  and  tested. 


Part  D:  Evolutions  of  fronts  in  flows  governed  by  nonlinear  equations  ("burning  surfaces") 

During  the  summer  of  1988  we  started  investigations  of  flows  generated  as  solutions  of  the  equation 

fofoy.O  &*(x,y,t)  &z(x,y,t)  r  32z(x,y,Q12  . 

3f  "  dx2  dy2  _L  dxdy  J 
(x,y)  e  £5,  t  €  (O.oo) 

z(x,y,0)  =  z„(x,y),  (x,y)  e  £2,  z(x,y,t)  |jo  =  <K*.y.O.  t  e  (0,<»), 
where /,  z0  and  $  are  given  functions. 

This  research  is  motivated  by  several  facts.  Fust  of  all,  under  some  assumptions  we  can  show  that  as 
t  -*  oo  the  solution  z(x,y,r)  converges  to  the  solution  of  the  steady  state  equation 


(6) 

(7) 
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32z  d2z 
dx2  dy2 


where  the  right  hand  side  /„  =  lim if  ( x,y,t )  as  t  — » ».  This  is  exactly  the  equation  (1).  Thus,  we  have  an 
alternative  method  for  constructing  solutions  of  (1),  (2). 


In  the  numerical  computations  this  approach  has  the  advantage  that  one  does  not  have  to  invert  the 
nonlinear  operator  M  but,  instead,  just  computes  its  value  on  the  current  update.  The  price  for  that  is  the 
need  to  do  a  lot  more  iterates  (in  time)  than  in  the  method  for  solving  (1),  (2)  as  described  in  part  B.  How¬ 
ever,  for  domains  with  complicated  geometry,  or  when  the  operator  on  the  right  hand  side  of  (6)  is  substan¬ 
tially  more  complicated  than  M,  this  may  be  an  acceptable  payoff  because  each  iteration  here  is  fairly  sim¬ 
ple  to  perform,  and  fewer  checks  along  the  way  are  needed. 

Another  factor  motivating  the  study  of  (6)  stems  from  the  observation  that  (6)  describes  a  flow  with  a 
front  moving  with  curvature  dependent  speed  (the  first  term  on  the  right  of  (6)  is  proportional  to  the  Gauss 
curvature  of  the  surface  front).  Similar  flows  were  used  in  physics  to  describe  combustion  models  (this  is 
the  origin  of  the  term  "burning  surface").  Flows  involving  Gauss  curvature  arise  in  several  applications; 
for  example,  the  rate  of  wear  of  certain  materials  is  proportional  to  the  Gauss  curvature. 

The  equation  (6)  is  also  interesting  in  itself  because  it  represents  a  flow  governed  by  a  genuinely 
nonlinear  differential  operator  in  spatial  derivatives. 

Many  of  the  ingredients  constructed  in  parts  B  and  C  were  utilized  in  developing  a  discretization 
scheme  for  (6).  We  have  written  a  code  for  numerical  solution  of  (6)  and  tested  it  on  a  number  of  exam¬ 
ples.  This  code  proved  to  be  very  useful  in  testing  various  conjectures  about  flows.  On  the  other  hand  the 
numerical  experiments  that  we  ran  suggested  a  number  of  interesting  properties  of  such  flows,  not  all  of 
which  have  been  known  even  for  classical  parabolic  equations.  The  results  of  this  investigation  are  in 
paper  "Evolution  of  nonparametic  surfaces  with  speed  depending  on  curvature,  I.  The  Gauss  curvature 
case".  A  special  case  of  such  a  flow  is  a  front  propagating  with  normal  speed  proportional  to  the  Gauss 
curvature  of  the  frontal  surface.  In  this  paper  sharp  estimates  of  the  asymptotic  rates  of  decay  are  esta¬ 
blished  and  these  results  are  applied  to  show  that  such  flows  "pick  up"  asymptotically  the  symmetries  of  the 
boundary  of  the  domain. 

Part  E:  Numerical  solutions  of  nonlinear  equations  and  computer  graphics. 

Since  the  iterates  in  numerical  approximations  of  solutions  of  (1),  (2)  are  represented  as  convex  hulls 
of  sets  of  points  in  /? 3 ,  it  is  natural  to  represent  those  iterates  as  graphs  of  polyhedral  surfaces  in  /?3.  Dur¬ 
ing  development  of  algorithms  for  numerical  solutions  of  PDE’s  it  is  extremely  useful  to  see  on  the  com¬ 
puter  screen  the  surfaces  corresponding  to  each  iterate  and  how  they  evolve  in  the  process  of  updating.  We 
could  not  find  appropriate  software  on  the  market  and  developed  our  own,  suitable  for  our  purposes.  The 
software  that  we  developed  generates  a  sequence  of  surfaces  on  the  screen.  Those  surfaces  are  colored 
depending  either  on  particular  physical  properties  of  the  solution  or  on  the  point  of  vision  of  the  observer. 

PartF:  Applications 

A  particular  nice  practica1  application  of  nonlinear  equations  of  Monge-Ampeie  type  is  'he  problem 
of  antenna  design.  In  this  problem  it  is  required  to  determine  a  reflecting  surface  such  that  for  a  given 
point-source  of  light  the  reflected  rays  cover  a  prescribed  region  co  of  the  far  sphere  and  the  density  of  the 
distribution  of  reflected  rays  is  prescribed  in  advance  as  a  function  of  the  reflected  direction.  It  is  assumed 
that  the  power  density  of  the  source  as  well  as  the  aperture  of  the  incidence  ray  cone  are  known,  and  the 
reflection  process  obeys  the  laws  of  geometric  optics.  In  this  form  the  problem  was  posed  by  Westcott  and 
others.  The  research  of  these  investigators  has  been  supported  by  Plessey  Radar,  Ltd.  for  many  years. 
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The  problem  admits  a  precise  mathematical  formulation  and  in  this  form  it  reduces  to  solving  the 
equation 

4P2det[vi;P  +  (P  -  [V-^+p2)eij] 

(1 VF  |2  +  P2)2det(ei;)  =fin  (8) 

nonlinear  boundary  condition  on  3co  (9) 

with  respect  to  the  unknown  function  P  (>0)  naturally  associated  with  the  problem;  here  (eiy)  is  the  matrix 
of  the  first  fundamental  form  e  of  the  unit  spheres  S2,  V  the  gradient  in  the  metric  e,  V1;- second  covariant 
derivatives  in  e,  and  /  the  prescribed  power  density.  The  condition  (9)  is  somewhat  complicated  to  be 
presented  here  without  considerable  expansion. 

With  the  support  of  AFOSR  we  started  our  investigation  of  the  problem  with  the  radially  symmetric 
(r.s.)  case,  that  is,  when  the  incidence  ray  cone  Cl  and  the  far  field  domain  to  are  circular,  the  prescribed 
density  of  reflected  rays  is  a  function  of  the  asimutal  angle  only,  and  the  reflecting  surface  is  sought  as  a 
surface  of  revolution.  Later  we  have  also  investigated  the  noradially  symmetric  case.  We  obtained  the  fol¬ 
lowing  results: 

•  the  r.s.  case  is  completely  resolved.  Simple  and  verifiable  necessary  and  sufficient  conditions  for 
solubility  of  (8)  and  (9)  are  given  and  it  is  shown  that  radially  symmetric  solutions  can  be  con¬ 
structed  explicitly  whenever  all  of  the  parameters  are  appropriate. 

•  it  is  shown  that  in  the  nonradially  symmetric  case  with  circular  incidence  ray  cone  and  far  field  and 
distribution  density  close  (in  certain  norm)  to  a  radially  symmetric  distribution  this  problem  admits 
two  classes  of  solutions  provided  the  data  satisfies  a  necessary  condition  expressing  the  energy  con¬ 
servation  law. 

We  also  began  investigating  the  following  inverse  problem.  Consider,  as  before,  a  reflector  antenna 
system,  consisting  of  a  point  light  source  O,  a  reflecting  surface  F  and  object  T  in  space  to  be  illuminated  in 
this  system.  Under  the  assumptions  of  geometric  optics  theory  it  is  desired  to  construct  the  surface  F  if  the 
positions  of  the  light  source  O  and  object  T  are  given,  and  the  power  distribution  is  a  prescribed  in  advance 
function  on  T.  In  addition,  the  aperture  of  the  incidence  ray  cone  is  a  design  parameter  and  it  is  also 
prescribed.  For  this  situation  we  developed  the  appropriate  equation  relating  the  phase  and  power  distribu¬ 
tion  generated  by  a  reflector  illuminating  the  target  In  our  considerations  the  target  is  allowed  to  be  an 
arbitrary  surface.  This  is  a  considerable  improvement  of  an  earlier  result  of  British  researchers  of  Brickell 
and  Westcott  Further,  we  have  given  general  conditions  for  a  given  function  to  be  the  phase  of  a  reflector 
illuminating  an  imbedded  surface.  Then  we  established  conditions  for  solvability  of  the  above  inverse 
problem  and  showed  that  under  natural  restrictions  one  can  construct  a  reflector  surface  which  will 
illuminate  a  given  surface  of  revolution,  coaxial  with  the  reflector,  with  prescribed  in  advance  varying  light 
intensity  and  aperture  of  the  incidence  ray  cone. 

n.  Singular  Boundary  Value  Problems 
The  nonlinear  boundary  value  problem 

+/<*,«)  =  0  (*) 
u  |r  =  $(x) 

on  a  region  fi  in  R*  with  boundary  T  is  a  problem  of  mathematical  interest.  The  special  case 
F  (x,u)  =  a(x)up  ,p  >  0,  is  known  as  the  Emden-Fowler  equation  and  is  of  considerable  interest  in  the 
applied  literature.  The  case  that  p  >  0  is  also  of  interest;  p  <  0  makes  the  problem  singualr  and  is 
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chaUenging  for  both  theoretical  and  numerical  reasons.  Many  of  the  computational  difficulties  are  already 
present  in  the  ordinary  differential  equation 

y”  +  a(x)yp  =  0 


with  the  usual  Sturm-Liouville  boundary  conditions.  Also  of  interest  is  the  problem  which  yields  radially 
symmetric  solutions  of  (*), 

y*  +  2-  +  a  ( x)yp  —  0 
y'(0)  =  0,y(l)  =  0. 

Early  in  the  grant  generalizations  of  both  the  Emden-Fowler  version  and  the  singular  version  were  investi¬ 
gated.  The  equation  took  the  form 

y"+f(x,y)=o 

with  general  Sturm  Liouville  boundary  conditions  and  with  conditions  on  /  that  included  either  the 
Emden-Fowler  equation  or  the  singular  equation.  (The  singular  equation  is  used  to  model  the  motion  of 
psuedo-plastic  fluids.  The  method  of  proof  involved  the  construction  of  iterative  schemes  and,  in  the 
singular  case,  a  fixed  point  theorem  for  decreasing  mappings  which  is,  in  its  own  right,  of  mathematical 
interest  For  the  non-singular  case  there  is  no  numerical  advantage  in  the  iterative  scheme  since  an  integra¬ 
tion  is  involved.  However,  in  the  singular  case  shooting  methods  are  not  appropriate  since  the  initial  value 
problem  may  not  be  well  posed.  The  singularity  effectively  prohibits  the  use  of  finite  differences  schemes 
as  well,  so  the  iterative  methods  take  on  new  importance.  Moreover,  iterative  schemes  lend  themselves  to 
parallel  computation.  The  system  was  programmed  and  the  iteration  technique  found  to  be  effective. 
Because  the  programming  is  so  simple  (no  large  ODE  solver  is  needed),  it  can  be  done  easily  on  very  sim¬ 
ple  personal  computers.  (Some  of  the  early  computations  were  done  on  an  Atari  to  illustrate  this  point.)  It 
is  clear  that  for  the  partial  differential  equation  larger  machines  will  be  needed.  Since  many  of  the  compu¬ 
tational  difficulties  are  present  in  the  boundary  value  problem  for  the  ordinary  differential  equation,  a 
parallel  version  has  been  implemented.  The  idea  is  to  use  one  processor  to  collect  the  data  and  assign  the 
task  (integrations,  comparisons,  output)  to  the  other  processors.  This  is  currently  in  progress  and  it  is  the 
hope  that  much  can  be  learned  which  will  be  useful  in  the  full  PDE  (singular)  equation.  The  code  for  a 
program  already  exists  and  u  has  been  obseieved  experimentally  that  the  speed  of  the  computations 
increases  almost  linearly  with  respect  to  the  number  of  processors  available;  the  computations  have  been 
carried  out  in  an  Encore  Mulitmax  computer  with  twelve  processors  but  the  program  will  run  with  any 
number  of  processors. 
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